Phase diagram and Debye mass in thermally reduced QCD 
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At temperatures well above the transition QCD admits a thermally reduced version 
in 3D, which reproduces the long distance physics. We analyze the phase diagram, 
point out the relevance of Z{3) symmetry in the location of the transition and 
suggest a way out to reconcile this with the data. Related to this symmetry is 
the existence of an observable, the Z(N) wall, or rather its 3D version, and discuss 
some of its advantages over other observables. 



1 Introduction 

One of the main challenges of thermal QCD is to geLreliable numbers. Though 
the gauge coupling may be small, Linde's argument tells us that perturbation 
theory will fail. The powerlike infrared divergencies one meets in perturbation 
theory will off-set the powers of the coupling constant. At what order in 
perturbation theory this will happen depends on the observable in question. 
For the free energy this happens when the static sector starts to dominate, 
and a simple dimensional argument shows this will happen at 0{g^). For the 
Debye mass Linde's phenomenon starts already at next to leading order. So 
the problem is certainly not academic! One should bear in mind that Linde's 
argument does not deny the existence of a perturbation series. It says that 
from a certain order on the coefficients are no longer obtained by evaluating a 
finite number of diagrams of a given loop order. 

So we are faced with evaluating non-perturbative effects from the three 
dimensional sector defined by the static configurations. It was realized some 
time ago □ that one could take the static part of the 4D action combined 
with induced effects by the non-static configurations. This theory gives at 
large distances the same physics as the 4D theory, and has the advantage 
of relatively straightforward lattice simulation^!. In section 2 we discuss the 
relation between the 4D and the 3D theory. In particular we show how the 
phase diagram of the 3D theory has a remarkable property: the curve of 4D 
physics, and the critical curve as determined by perturbation theory do coincide 
to one and two loop order. However, perturbation theory has no reason to be 
trustworthy in determining the critical curve, and this is probably the reason 
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why the fit to the numerical determination is problematic. 

In section 3 we discuss the physics of the domain wall in some detail. 

2 Effective 3D action and symmetry in 4D 

Construction of the effective action proceeds along familiar lines. In the case 
of QCD with Uf quarks its form is given by integrating out the heavy modes 
of 0(r): 

S3D = SyM, n=0 + Sind (1) 

The first term is the static sector of the pure Yang-Mills theory in 4D with 
coupling constant (73 = gVT. 

The second term in eq. Qmust contain the symmetries of the original QCD 
action, as long as they are respected by the reduction process. 

So we expect the induced action to be of the form: 

Sind — ^(^0) + terms involving derivatives (2) 

V{Ao) should be invariant under static gauge transformations, C, CP 
{Aq — ^o') ^^'^ t'^is reduces it to a sum of traces of even powers of Aq: 

V{Ao) = m^TrAl + \i{TrAlf + MTrA^^ + .... (3) 

Only one independent quartic coupling survives for SU(2) and SU(3). We 
take it to be TrA^. Note that we lost a symmetry present in the 4D action for 
gluons alone, and less and less conserved when quarks get lighter and lighter: 
Z(N) symmetry. 

Remember from the lattice formulation of pure Yang-Mills that one can 
multiply at a given time slice in the original 4D action all links in the time 
direction with a factor exp±i This will not change the form of the action, 
but will change by the same factor the value of the Wilson line P wrapping 
around the periodic time direction: 

P(Ao) =7'expi j AodT (4) 

Clearly in eq. ^ this symmetry has gone. Apparently the reduction process 
does not respect Z{?>) symmetry! The reason for this is twofold: 

i) the reduction process does not include the static modes. 

ii) the values of Aq/T in the effective action are order g, whereas the Z(N) 
symmetry equates the free energy in Aq/T and Aq/T + 0(27r/3). 

To understand this better - and to prepare the way for the discussion of 
the domain wall observable in the last section |^ - we recall some familiar facts 
in 4D for SU(3). 
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2.1 Z(3) symmetry and domain walls in 4D gauge theory 

The free energy J7 as a function of the Wilson line invariants TrP, TrP^ is 
naturally defined through: 

e^p-YEi^lllA ^ f DAoDAS{ti-T7P)d{t2~TW^) exp-^^ (5) 
T J g-^ 

where TrP is the normalized space average of the trace over the volume V . A 
natural parametrization of the parameters ii and t2 suggests itself: define the 
phase matrix expiC with C being a traceless diagonal 3x3 matrix with entries 
G, , (i = 1, 2, 3) and G, = 0, because we have SU(N), not U(N). 

Consider pure Yang-Mills. A gauge transformation that is periodic modulo 
a phase in Z(3) will only change the arguments in the delta functions in eq. |[ 
Hence the potential U has degenerate minima in all points of the C-plane, 
where expiC = 1, or exp±i27r/3. This is called Z(3) symmetry (and the 
degeneracy is lifted by the presence of quarks). 

This statement is independent of perturbation theory. In fact the potential 
in eq. |^ has been computed in perturbation theory including two loop order. 
And this potential includes the static modes. Propagators acquire a mass 
proportional to the phases C, because it acts like a VEV of the adjoint Higgs 

Hence, for small C, eventually Linde's argument will apply and the per- 
turbative evaluation becomes impossible. 

For SU(3) the direction in which the Wilson line phase causes minimal 
breaking is in the hypercharge direction G — ^diag{q, q, —2q). Minimal break- 
ing means the maximal number of unbroken massless excitations, that do not 
contribute to the potential. Hence this is at the same time the valley through 
which the system tunnels from one minimum to the next. In this "q- valley" 
the combined 1 and 2 loop result is exceedingly simple: 

[/(I) + C/(^) = ^-^T\N - 1) (l - 50) ,^(1 - qr (6) 

For use in the reduced theory we isolate the static part of the one and two 
loop contribution in the q- valley from eq. |6[ 

(U^'^ + U^'^) = -T\N -1)^4- f^'?' + ^7^2^ (7) 
V / (n=o) 3 \ [AttY J 

Note that the two loop contribution is quadratic in q in contrast to the one 
loop which is cubic. The two- loop cubic part in eq. ^ comes from a combination 
of static and non-static modes. 
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If we prepare the 4D system conveniently this symmetry will give rise to 
domain walls. Profile ancLenergy of these wall have been computed semi- 
classically a long time agco. The method of twisted boundary conditions 
triggers walls and is most economic computerwise. We will discuss them in the 
context of the lattice formulation in section |3[ Be it enough to mention that 
these boundary conditions force the Wilson lines to change by a Z(N) phase in 
going from one side to another side of the volume in some a priori fixed space 
direction. This will trigger a wall profile for the loop in this direction. 

It is the long range behaviour of this profile that contains the information 
on the Debye mass. To one loop order this behaviour comes entirely from the 
slope of the potential, see above. But to two loop order we have to take the 
one-loop renormalization of the gradient part of the Wilson line phase into 
account, and this suffers the Linde effect: there is an infinity of many-loop 
diagrams contributing to the gradient part. So to next to leading order there 
are already non perturbative effects in the long range tail of the wall, and hence 
in the Debye mass, as we mentioned earlier. 

On the other hand we know that the effective 3D action correctly repro- 
duces the large distance behaviour of the 4D theory. So a 3D projection of the 
twist should produce a wall with the same tail as the 4D one. The inside of 
the wall in both formulations may be quite different but the inside is anyway 
computable by perturbation theory. 



2.2 3D action and 4D physics 

The parameters of the 3D theory (m^ and A = Ai + A2 for SU(3)) in eq. || can 
be calculated in perturbation theory by integrating out all modes in a path 
integral except the mode A^(a;, n = 0). To one loop order we have the well 
known result for the Debye mass and for the four point coupling A. All higher 
order terms have a coefficient zeroB. To two loop order one has to take care not 
only of the two loop graphs, but also of the 1-loop renormalization of the three 
dimensional gauge coupling 33 and the renormalization of the Aq field in the 
gradient terms. The latter renormalization is taking care of gauge dependence 
in the two loon graphs. 

The results in the MS scheme is that both parameters are expressed in the 
renormalized 4D coupling g{fi) where /i is the subtraction point. Eliminating 

the 4D coupling gives for the dimensionless quantities x = ^ and y = ^ the 

33 93 

result for N=3: 

3 3 

^2/4D = + ^x) (8) 
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whereas for N=:2: 

2 9 

xy4D = 7-^(1 + ^x) (9) 

Note the absence of exphcit /i dependence in this relation. The variable x 
has a ^ dependence such that as T becomes large x becomes small. 

In conclusion, it is along this line that we have to simulate the 3D system, 
in order to get information about the 4D theory. Before we do this, we still 
have to settle an important question: where are - in the xy versus x diagram 
- possible phase transitions? 



2.S Phase diagram of the 3D theory 

To get the phase diagram we must first decide what order parameters to take. 
In the case of SU(3) there are two: TrA^ and TrA^. Strictly speaking, only 
the latter is an order parameter, since it flips sign under C. We will study the 
analogue of eq]|: 

e^p-VSeff{D,E)= J DAS (glD - TrA^j 5 [qIE ~ TrA^j exp ~5 (10) 

Again as for the Wilson line we parametrize D and E in terms oi D = Tr [C^] 
and E = Tr [C^] respectively. Let us first state the result one gets for Seff to 
one and two loop order: 

t/(n = 0) , , , 

'-'eff — Tj; one and two loop only (11) 

The one and two loop result equals the static part of the 4D Z(3) poten- 
tial,eq. |^! This static part was explicitely written in the q- valley, eq.0. It has to 
be added to the tree result and one gets in terms of the dimensionless variables 
X and y for N= 2 or 3 colours, absorbing a factor 27r in q: 

The question is now: for what values of x and y we have degenerate minima 
for q? Keeping only the 1 loop result cubic in q we see that it must be of 
the order of magnitude of the quartic term of the tree result to get a second 
degenerate minimum. So q must be of O(^) in that minimum. Thus the 
quadratic two loop result contributes 0{x) less. 

From eq.|l^ we find the potential develops two degenerate minima for N=3 
when: 

3 3 

a;2/c = ^(1 + 2^;) (13) 
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For N=2: 

^?/^ = ^(l + ^^) (14) 

This is important: slope and intercept of the physics hne ^ are identical 
with those of the critical line ^ at least if we can take the low orckiLloop 
results for the critical line seriously. This was numerically found in reft3El The 
intercept equality is just due to the Z(N) potential in 4D and the effective 
potential Seff in 3D being identical to one loop. But to two loop order this 
simple explanation is no longer true. The cubic term in eqj^ is appearing also 
in the two loop result, but not in the two loop result for the 3D effective action. 
It is however true that also in 2 loops the leading contribution is the static part 
of the Z(N) potential, eqj^. 

2.4 Saddle point of the effective potential in 3D 

In this subsection we will investigate in more detail the computation of the 
3D effective potential. The saddle point is found by admitting Aq fluctuates 
around a diagonal and constant background B: 

Ao=B + Qo (15) 

whereas the spatial gauge fields fluctuate around zero: 

A, = g,. (16) 

One then goes through the usual procedure of expanding the effective action 
|To| . The equations of motion fix the background B to be equal to the matrix C, 
and the part quadratic in the fluctuations will not contain any reference to the 
Higgs potential V{Aq). This is clear because the quadratic constraint tells the 
mass term not to fluctuate. Only the Higgs component parallel to C , TrCQo, 
has a mass term due to the Higgs potential, AXTrC^. So apart from this the 
quadratic part comes entirely from the static part of the 4D action. We can 
make a convenient gauge choice, namely the static form of the covariant gauge 
fixing: 

Sgf = Tr {[tgsB, Qo] + dkQkf (17) 

This gives propagators which are precisely the static version of the prop- 
agators appearing in the Wilson line potential |5| Only the component Qo 
parallel to C is the exception: its propagator has a mass from the Higgs poten- 
tial and can be written as the the sum of the static propagator and a remaining 
part ( "massive" ) containing the mass term; 

111 

(18) 



p2 p2+AXTrC^ p2 
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The static propagator dominates in diagrams over the rest. The massive 
propagator will give rise to half integer powers of x in the perturbative expan- 
sion of the potential; gauge couplings contribute 0(1) in dimensionless units, 
whereas Higgs couplings contribute 0{x). 

As long as we are interested in intercept and slope of the critical curve, it 
follows that only the static part of the Feynman rules contributes. 

Hence the result ^ 

Let's from now on work in the q- valley where we evaluate the effective 
action |l2|. 

Then two remarks are crucial: 

i) The broken minimum occurs for q = 0{l/x) . Power counting then reveals 
that from 0{x^^'^) on an infinite number of diagrams contributes to each order. 

ii) From five loop order on, the potential starts to develop poles in q = 0. 
We are bringing this up, because insisting on the low order result |l^ and 

fitting numerically the coefficients of x^/^ and higher order giveS|-an unexpected 
result: the numerical coefficients are orders of magnitude largefl than the first 
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f3. In fig. 1, taken from ref. u, the situation is shown. Only for very 
small X the critical and the 4D physics line are allowed to become tangent. It 
seems that this constraint affects the quality of the fit. Dropping it altogether 
necessitates numerical determination of transition points at x < 0.04. 
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Figure 1: Phase diagram in the SU(2) case, from ref. 9. The straight line is the 4d 3d 
curve of eq.^. The thick line is a 4th order fit to the data. The dashed line marks the region 
where the transition turns into a cross-over. 
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3 Debye mass from a 3D domain wall 

After this long discussion of where the physics hnc hcs with respect to the 
critical curve we have to come to grips with the domain wall method. |— . 

The idea here is extremely sirrmle and has been explained elsewhercO. 
Twisted boundary conditions in 4Dll3 have a very simple and intuitive form 
in the reduced theory. Remember that a twisted plaquette in the time-space 
direction is of the form Tr(l — n,U{P)), with = expi2TT/N. 

Thus intuitively one would say that all one has to do in the reduced action 
is to modify the kinetic part of the Higgs field by the twist, because that's 
what the plaquette in the time space direction is reducing to. 

In the next subsection we work out this idea in more detail. 



3. 1 Construction of the wall 

In this section we want to make rniire precise the action that defines the wall. 

We follow the notation of refE2l, specifically that of hep-lat/9811004 and 
write the kinetic part of the action as: 



X X,j 



(19) 



where a; is a vector with three components {x,y,z). 
Consider the following expression: 



1 



1 



-Ue Tr 



^iaA(x,yfi)jj^-iaA(x,yS)jj+ 



If aA is small we get: 
1 



X = -a^Tr 
N 



\a\0) + \a\1)~ A{0)UA{l)U+ 



(20) 



This is precisely the kind of expression that appears in eq. (19). /^From this 
follows the expression for the modified kinetic energy in the plane (a;,y,0): 



\ (A^x) + A^ix + ae,)) - A(f)L/,(f)A(f + ae,)f/;(f) 
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(21) 



So all we need is to put a twist O G ^(N) in order to get a wall: 



pwall 



+ 1 1 - M^K" [e^"^("'^'°)«73(a;,t/,0)e--^(-'^'i)C/3+(x,y,0) 
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(^,j)^(0,3) 



{A^{x) + A'^{x + aej)) - A{x)Uj{x)A{x + aej)U^ (x) 



(22) 



What is now the actual value of a to use? We recover the kinetic term in 
the continuum if wc relate the field A on the lattice to the field Acont in the 
continuum by the relation: 

Acont 



A = 
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This is not the usual normalization for the lattice fields. Usually we have: 
Aiatt = ags-A-cont and Aiatt ^ in the continuum limit. 

Here this is not anymore the case. Remember that to expand the modified 
action wc had to suppose that aA was small. To enforce this condition it seems 
natural to put : a = ag^; in this manner terms of the kind e'"^ become e'^^a"^. 
That is to say, they become of the usual sort : e'"^=*^=°"' . 

With this choice the term in the exponential indeed goes to zero as the 
lattice spacing goes to zero, so: 

_ 2 6 
a = ag^ = - 

In the end we obtain as final expression for the kinetic part of the action 
supporting the wall: 



r*wall 



{A^{x) + A^{x + aej)) - A{x)Uj{x)A{x + aej)Uf{x) 



+/3 



T 11 



VL Tr 



J§Ai^,y,o)u^^^^ 0)e~'t^(">«'i)?73+(x, y, 0) 



(23) 
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3.2 Excitations of the wall 



Now the system with the wall is defined by adding the 3D gauge field action 
and the Higgs potential V(A) to eq. |2^. Let us call the resulting twisted action 
St. 

Both twisted and untwisted action have periodic boundary conditions. 
When we compute the average of an observable O in the twisted box we aver- 
age the observable over the {x,y) plane at the point z, written as 0{z), and 
compute in the twisted box (action St). It is quite trivial to relate this average 
to the correlation of the wall and O in the untwisted box (action S): 

(0(i)>s. = (exp -(5t - S)0(r))s (24) 

There is no difference between the two actions except at z = 0, at the location 
of the wall. 

The twist is C and P odd, but T even. This means we can expect a sigiial 
for the Debye mass by taking any observable O C odd (a necessary conditioiH). 
Whatever operator gives the lowest mass in the correlation |2j is the preferred 
one. Thus one and the same updating with the twisted box can be used for 
various operators. 



4 Conclusions 

Once we know the 4D physics line we can do a simulation of the twisted box 
with some convenient observable, and measure the mass through eq. Care 
should be taken, as emphasized by Kajantie et al. that we start in the 
symmetric phase and then move to the 4D physics line. In so doing we will 
stay on the physical branch of the hysteresis curve for the mass, that we will 
meet when crossing the transition curve. 

Nethertheless our discussion of the location of the critical curve underlines 
the importance to know wether the 4D physics line lies for small x in the 
symmetric phase or in the broken phase. 
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